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ABSTRACT 


The  Airborne  Laser  Laboratory  (ALL)  is  being  developed  by  the 
Air  Force  to  investigate  the  integration  and  operation  of  high 
energy  laser  components  in  dynamic  airborne  environment  and  to 
study  the  propagation  of  laser  light  from  an  airborne  vehicle  to 
an  airborne  target.  The  ALL  is  composed  of  several  systems;  among 
these  are  the  Airborne  Pointing  and  Tracking  System  (APT)  and  the 
Automatic  Alignment  System  (AAS).  This  report  presents  the  results 
of  performing  a  time  domain  dynamic  simulation  for  an  integrated 
beam  control  system  composed  of  the  APT  and  AAS. 

The  simulation  is  performed  on  a  digital  computer  using  the 
MIMIC  language.  It  includes  models  of  the  dynamics  of  the  system 
and  of  disturbances.  Also  presented  in  the  report  are  the  rationales 
and  developments  of  these  models. 

The  data  from  the  simulation  code  is  summarized  by  several  plots. 
In  addition  results  from  massaging  the  data  with  waveform  analysis 
packages  are  presented.  The  results  are  discussed  and  conclusions 
are  drawn* 
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I.  INTRODUCTION 


The  Airborne  Laser  Laboratory  is  composed  of  several  systems. 
Among  these  are  the  Airborne  Pointing  and  Tracking  System  (APT)  and 
the  Automatic  Alignment  System  (AAS).  The  primary  functions  of  the 
APT  are  (1)  performance  of  precision  tracking  and  (2)  stabilisation 
of  the  optical  line  of  sight  (LOS)  to  the  target.  The  primary 
function  of  the  AAS  is  to  compensate  for  static  misalignments  of 
the  laser  beam;  a  secondary  function  is  to  aid  the  APT  in  the 
stability  of  the  LOS.  Two  causes  of  instability  of  the  LOS  are 
base  motion  disturbance  (BMD)  and  target  motion.  In  order  to 
reduce  the  effects  of  BMD  on  the  stability  of  the  beam  LOS  studies 
were  initiated  to  determine  the  feasibility  of  using  the  Automatic 
Alignment  System  in  conjunction  with  the  APT,  for  beam  control. 

(Ref.  1,  2). 

The  scheme  devised  for  reducing  BMD  involves  feeding  the  stab 
gyro  output  through  a  100  Hz  low  pass  filter;  then  to  the  AAS  and  to 
the  image  motion  compensation  (IMC)  input  of  the  tracker.  (See 
Figure  1).  The  100  Hz  low  pass  filter  minimizes  errors  due  to  high 
frequency  gyro  resonances.  Feeding  the  gyro  output  through  the 
100  Hz  filter  to  the  image  motion  compensation  and  to  the  AAS  pro¬ 
duces  a  transf-'r  function  between  the  BMD  input,  6^^  »  anc* 
beam  angle,  6  ,  that  is  approximately  zero  for  frequencies  below 

D 

100  Hz.  For  frequencies  significantly  greater  than  100  Hz  the 
affect  on  by  0_._.  is  practically  unaltered  by  the  feedforward 

D  oMD 


I 


occurs  in  the 


scheme.  Fortunately  most  of  the  energy  of 

“  BMU 

frequency  range  below  100  Hz.  Furthermore,  the  pure  integration 
between  9^^  and  will  greatly  reduce  the  effects  on  the 

beam  angle  of  signal  components  in  the  greater  than  100  Hz  range. 

The  above  scheme  is  shown  in  [2]  to  greatly  improve  BMD 
rejection  but  does  not  improve  the  capability  of  the  system  to  track 
targets  with  significant  motion  relative  to  the  APT/AAS.  In  order 
to  provide  this  capability  it  was  decided  to  feed  the  tracker  error 
to  the  input  of  the  AAS.  This  results  in  an  overall  increase  in 
the  bandwidth  of  the  system  by  a  factor  of  twenty.  The  major 
disadvantage  of  this  scheme  is  that  the  effect  on  0_  due  to 

D 

tracker  noise  is  amplified. 

The  tracker  noise  and  the  target  signal  0g  are  modeled  as 
entering  at  the  same  point.  (See  Figure  1).  Thus,  the  combined 
beam  control  system  (APT/AAS)  has  the  same  filtering  affect  on  both 
signals.  Fortunately  the  power  density  spectrums  (PSD's)  of  0g 
and  0T  differ.  Typically  for  a  "flyby"  target  the  PSD  of  0g  is 
low  frequency  dominated  where  the  PSD  for  0^  is  approximately 
flat  up  to  30  Hz.  The  transfer  function  P(s)  in  Figure  1  is 
inserted  to  provide  compromise  filtering  between  0^  and  0^  so  as 
to  reduce  the  beam  error,  (0  -  0  ). 

b  O 

In  an  ideal  system  the  transfer  function  between  9  and  0 

D  5 

should  be  unity.  However,  because  of  the  system  dynamics  a  unity 
transfer  function  is  impossible.  The  beam  angle  is  a  function  of 


0 


S 


9_  ,  and  0_,_,  »  consequently,  the  beam  error  is  also  a 

l  15  MU 


function  of  these.  For  the  system  to  perform  high  precision 
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tracking,  the  beam  error  must  be  small.  In  the  frequency  domain 
this  corresponds  to  requiring  the  PSD  of  the  beam  error  to  be  as 
small  as  possible  over  the  frequency  spectrum. 

The  design  criterion  for  the  filter  then  is  the  P(s)  that 
minimizes  the  root  mean  square  (RMS)  of  the  beam  error.  In  [3]  a 
computer  code  for  performing  this  minimization  is  presented.  This 
code  was  used  to  design  optimal  P(s)'s  for  several  cases  of  "flyby”  type 
scenarios.  (See  Figure  2).  The  results  are  summarized  in  Table  I. 
Data  is  presented  for  comparing  the  no  feed  forward  to  the  feed 
forward  with  P=1  case  and  the  optimal  P's  cases  and  for  comparing 
the  P**l  to  the  optimal  P's  cases.  The  P=1  case  is  presented  since 
this  is  the  logical  choice  for  P  without  optimization.  In  addition 
improvements  of  the  optimal  P's  to  the  P*1  cases  are  presented. 

This  table  shows  that  simple  optimal  lead-lag  type  filters  in 
the  feed  forward  scheme  can  produce  marked  reductions  in  the  RMS 
beam  errors. 

The  results  obtained  in  [3]  show  the  feasibility  of  controlling 
the  beam  angle  with  a  high  degree  of  accuracy.  However  the  figures 
of  merit  are  based  on  RMS  values  of  error.  , RMS  errors  do  not  provide 
any  information  concerning  the  dynamics  of  the  errors,  e.g. , 
maximum  and  minimum  values  and  the  points  in  time  of  occurance. 

In  addition,  in  the  study  in  [3]  the  sampling  effects  of  the  tracker 
were  Ignored.  Thus,  in  order  to  validate  the  results  obtained  in 
[3],  a  time  domain  simulation  of  the  feed  forward  scheme  is 
presented  in  this  report. 


The  simulation  is  implemented  in  the  MIMIC  Language.  (Ref.  4). 
The  development  process  for  the  simulation  involves  developing 
models  of  the  signals  and  the  blocks  of  Figure  1.  These  developments 
are  presented  in  the  two  sections  that  follow.  Following  these  is 
the  results  from  executing  the  simulation  code,  along  with 
analyses  of  the  results  and  conclusions. 
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II.  SIGNAL  MODELS 


Figure  1  shows  three  sources  of  external  signals  entering  the 

APT/AAS  pointing  system.  These  are  the  target  signal  0_  ,  the 

s 

tracker  error  0_  and  the  base  motion  disturbance  0BWT.  .  In 

i  JDflU 

[3]  PSD's  of  each  of  these  were  used;  however,  here  time  domain 
representations  of  these  signals  are  needed.  These  time  domain 
models  are  developed  in  the  following. 

A.  TARGET  SIGNAL  MODEL 

A  "flyby"  type  scenario  as  shown  in  Figure  2  is  assumed  for 
the  target  signal.  The  equation  that  describe  the  angle  of  the 
target  w.r.t,  the  APT/AAS  is 

0-  -  tan"1  (1) 

S  u 


where 


X  ■  Xq  +  vt  . 


(2) 


In  (1)  and  (2)  Yq  is  the  minimum  approach  of  the  target  to  the 
APT/AAS  pointing  system,  v  is  velocity  (assumed  to  be  constant) , 
t  is  time,  and  Xq  is  the  initial  x  value  where  tracking  is  to 
be  initiated.  Tracking  is  assumed  to  begin  two  seconds  before  the 
target  reaches  its  closest  approach  to  the  APT/AAS  and  continues 
for  four  seconds  afterward.  Therefore,  the  tracking  duration  is 
always  four  seconds,  but  the  point  where  tracking  starts  depends 
upon  v  and  Yq  . 


In  order  to  minimize  transient  effects  due  to  the  initiation 
of  tracking  it  is  assumed  that  an  automatic  handoff  system  will 
exist,  i.e,  the  initiation  of  tracking  will  be  smooth.  Without 
the  automatic  handoff  the  derivative  of  the  target  signal  appears  «>s 
typified  in  Figure  3(a).  At  the  times  where  tracking  starts  Tq 
and  ends  the  derivative  changes  discontinuously.  The  discon¬ 
tinuity  at  Tq  causes  transient  effects  that  can  be  detrimental 
to  good  tracking  performance.  When  tracking  starts  it  is  best 
that  the  beam  angle  be  positioned,  moving  and  accelerating  the 
same  as  the  LOS  to  the  target.  An  automatic  handoff  system  will 
provide  these  initial  conditions.  In  order  to  simulate  the 
automatic  handoff  system  a  fourth  order  polynomial  spline  is 
attached  to  the  target  signal  (thus,  generating  a  fictitious  target 
signal  for  time  less  than  TQ  ).  The  requirements  on  the  spline 
are 

0'(To-4)  -  0 
6' (Tq-4)  =  0 
0'(TO)  =  0(4) 

0'(TO)  -  0' (4) 

O' (Tq)  =  0(4) 

For  simplicity  in  the  following  it  is  assumed  that  Tq-4  *0.  In 
essense  the  target  signal  will  be  "faked"  by  the  spline  for 
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(a)  Typical  Signal  Rate 
e(T) 


(b)  Target  Signal  Rate  with 
Spline  on  Left 


Figure  l\  Typical  Target  Signal 
Angular  Velocity 
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0<t^4  seconds  and  from  4^t</  tracking  of  the  target  will 
occur*.'  The  Impact  on  the  derivative  of  the  target  signal  by  the 
spline  is  illustrated  in  Figure  3(b). 

B.  TRACKER  NOISE  MODEL 

In  [2]  the  tracker  noise  is  approximated  as  having  a  flat  PSD 

up  to  30  Ha  and  having  an  RMS  value  of  3.1  rad.  In  actuality  the 

tracker  noise  is  a  train  of  uniformly  occurring  random  pulsea,  e.g., 

see  Figure  4.  The  amplitude  of  each  pulse  la  assumed  to  be 

2 

normally  distributed  with  aero  mean  and  a  variance  of  a  .  The 
power  spectral  density  of  a  waveform  of  this  type  is 

♦  (u)  ■  o2T  (3) 

<u>t/2  r 

where  T  (sampling  period)  is  the  duration  of  each  pulse.  For  the 

2  -12 

tracker  being  considered  in  this  study  a  •  9.61x10  and 
T  0.01666s. 

A  plot  of  the  power  spectrum  density  of  the  model  used  in  this 
study  is  shown  in  Figure  5.  The  PSD  is  approximately  flat  up  to 
30  Hr  as  specified  above.  Thus,  for  generating  this  signal  in  the 
MIMIC  coda  all  that  is  required  is  to  sample  and  hold  every  T 

seconds  a  normally  distributed  waveform  with  zero  mean  and  a  variance 

,  2 
of  o 


*  The  discontinuity  at  the  end  of  tracking  T_  is  not 
considered  detrimental  since  it  occurs  after  tracking  has  ceased. 
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C.  BASE  MOTION  DISTURBANCE  MODEL 


In  [1),  [2],  and  [3]  an  experimental  PSD  for  BKD  was  utilised. 

In  order  to  generate  a  tine  domain  signal  with  the  characteristics 

of  the  experimental  PSD,  a  normally  distributed  train  of  pulses 

2 

with  xero  mean  and  variance  of  c  is  generated  and  filtered  to 
obtain  the  desired  PSD  characteristics.  The  PSD  of  such  a  waveform 
is 

*BMD<U>)  "  o2T  ~ n-(--T-^  |G(ju)|2  (A) 


where  T  is  the  pulse  durations  and  G(jw)  is  the  filter.  By 
trial  and  error  the  filter  that  produces  a  PSD  approximately  equal 
to  that  of  the  previous  studies  is 


G(ju>) 


iS  +  l 
9  X 


fjSL +J 
(  200  ^ 


(5) 


•.-A 


when  T  «  0.00628319  and  o  *  l.A32Axl0  "  .  A  plot  of  the  PSD 
is  shown  in  Figure  6.  Circles  on  Figure  6  indicate  points  from 
the  experimental  PSD  used  in  [3].  A  close  correlation  between 
the  data  is  seen. 
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Figure  5,  PSD's  of  Tracker  Noise  Figure  6.  PSD's  of  BMD 


III.  DEVELOPMENT  OF  THE  MIMIC  SIMULATION 


In  order  to  expedite  the  development  of  the  MIMIC  simulation. 
Figure  1  is  redrawn  in  a  more  genet el  form  in  Figure  7.  On  this 
figure  MIMIC  variable  assignments  have  been  watte-  e.g. ,  tracker 
error  Is  denoted  as  TRKERR  and  base  motion  disturbance  is  BMD, 

The  transfer  functions  used  in  the  development  of  the  MIMIC 
simulation  are  as  follows: 


G1(s) 

Vs) 

G3(s) 

Ga(s) 


e-0.003s 

33.02s+350 

s 

1 

s 

_ 265.3s  +  10A _ 

1. 274xl0~8s4+3.568xl0”6s3+4. 715xl0~3s2+s 


G5(s) 

G6(s) 


1 

s 

1987s-H.248xlQ6 

2 

s 


H^s) 


1.0 _ 

1.592xlO”3s-tl.O 


(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 


P(8)  depends  on  the  case  being  considered,  i.e,  cither  zero  for  no 
feed  forward,  unity,  or  anyone  of  the  optimal  P(s)'s  from  Table  I. 
In  addition  for  the  no  feed  forward  and  stab  not  fed  forward  cases 
the  100  Hz  filter  is  set  to  zero. 


. s», r  , 


It  should  be  noted  that  G^(s)  Is  different  from  that  used  in 
[3j.  In  this  previous  study  G^(s)  is  modeled  as  a  three  milli¬ 
second  delay  in  cascade  with  the  transfer  function  of  a  zero  order 
hold.  This  is  an  approximation.  In  actuality  the  tracker  involves 
a  sampling  process.  In  order  for  the  MIMIC  simulation  to  be  a 
more  accurate  representation  of  the  actual  system  the  equivalence 
in  this  report  of  G^(s)  is  a  3  millisecond  delay,  whose  transfer 
function  is  (6),  in  cascade  with  a  sample  and  hold,  i.e.  the  cascade 
combination  of  G^(s)  and  the  block  labeled  S/H  in  Figure  7. 

Using  Figure  7,  converting  each  transfer  function  to  a 
corresponding  set  of  differential/algebraic  equations,  and  using 
the  signal  models  of  the  preceding  sections  the  MIMIC  code  presented 
in  the  Appendix  was  developed  for  simulating  the  APT/AAS  beam 
control  system.  The  code  has  been  commented  for  ease  in  following. 
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IV.  SIMULATION  RESULTS 

The  code  used  in  Appendix  A  is  used  to  simulate  the  APT/AAS 
beam  control  system  for  four  c.ase3  of  the  flyby  type  scenario,  in 
particular,  Yq  m  500m,  1000m,  2000m,  and  5000m.  In  addition  five 
different  beam  control  configurations  are  simulated:  (1)  no 
feed  forward  (2)  tracker  error  not  fed  forward,  (3)  stab  error  not 
fed  forward,  (4)  complete  feed  forward  with  P(s)  *  1,  (5)  and 
complete  feed  forward  with  the  optimal  P(s).  In  the  second  con¬ 
figuration  the  stab  error  is  still  fed  to  the  AAS  and  fed  back 
to  the  IMC.  In  the  third  configuration  the  path  from  the  stab 
error  to  the  AAS  and  IMC  are  severed  and  the  optimal  P(s)  is 
used  in  the  tracker  error  to  AAS  path. 

Target  signals  of  the  four  cases  of  the  flyby  type  scenario 
are  shown  in  Figure  8.  The  plots  include  the  spline  signals  (0  to 
4  seconds)  plus  the  assumed  target  signals  (4  to  8  secondsl).  As 
expected  the  closer  the  target  comes  to  the  APT/AAS  the  greater 
the  angular  velocity  and  acceleration  of  the  target  with  respect 
to  the  APT/AAS. 

Results  from  the  simulation  for  determining  the  performance 
of  the  five  configurations  and  for  the  four  cases  of  the  flyby 
scenario  are  shown  in  Figure  9-13.  Tracker  noise  and  base  motion 
disturbances  are  included  in  all  simulations.  Each  figure  shows 
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Figure  13.  Beam  Error  for  Optimal  Filter  Configuration 


the  bean  error  in  v>rad  versus  time  in  seconds  for  the  3  to  8  second 
range.  This  time  range  is  seiected  because  the  main  interest  is 
the  tracking  ability  of  the  APT/AAS  of  the  actual  target  signal, 
not  of  the  spline. 

In  Figures  9  and  IQ,  large  beam  angle  errors  (greater  than 
±50  yrad)  occur  for  the  500m,  1000m  and  2000m  cases.  This  is 
caused  by  the  target  having  significant  relative  motion  to  the 
APT/AAS  (both  velocity  and  acceleration).  For  the  500m  and  1000m 
cases,  the  errors  drastically  Increase  before  and  after  the  time 
of  closest  approach  by  the  target  (t“6s) .  This  is  attributed  to 
the  target  angle  having  large  angular  accelerations  before  and 
after  t*6s.  On  the  other  hand,  the  2000m  and  5000m  cases  show 
the  errors  tending  to  montonically  decrease.  Contrary  to  the 
500m  and  1000m  cases  these  also  show  the  randomness  due  to  BMD  and 
tracker  noise.  Of  course  this  is  because  the  first  two  cases  are 
dominated  by  target  motion  errors,  which  are  many  times  greater 
than  errors  due  to  the  tracker  or  BMD. 

It  should  also  be  noted  that  for  the  5000m  case  the  error 
abruptly  decreases  just  prior  to  t*4s.  This  is  due  to  the  spline 
generating  a  large  fictitious  acceleration  prior  to  the  initiation 
of  tracking.  In  fact  for  the  5000m  case  the  APT/AAS  exhibits  fair 
tracking  capabilities  for  the  no  feed  forward  and  tracker  error 
not  feed  forward  configurations,  in  particular,  the  errors  range 
between  ±20  yrad.  For  the  other  cases  of  these  configurations  the 
errors  are,  in  general,  too  large. 
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In  Figures  11-13,  the  results  from  simulating  the  different 
cases  of  the  stab  error  not  feed  forward,  unity  feed  forward  and 


optimal  filter  feed  forward  configurations  are  summarised.  In 
comparing  the  corresponding  cases  of  any  of  these  configurations 
to  those  of  Figures  9-10,  it  is  easily  seen  that  the  tracking 
capabilities  have  been  definitely  augmented.  Of  course  the 
commonality  of  these  three  configurations  is  that  all  feed  the 
tracker  error  signal  to  the  AAS  in  some  form. 

Both  the  optimal  filter  and  stab  error  not  fed  forward  cases 
feed  the  tracker  error  signal  to  the  AAS  through  the  optimal 
filters  designed  in  [3].  Both  of  these  cases  provide  less  beam 
error  in  all  four  cases  of  the  flyby  type  scenario  than  the 
unity  filter  case.  This  was  expected  since  the  optimal  filter 
tends  to  compensate  for  the  tracker  not  being  perfect  (continuous 
output  with  a  unity  transfer  function) . 

The  stab  error  not  fed  forward  and  the  optimal  filter  config¬ 
urations  are  about  equal  in  performances  for  the  500m  and  1000m 
cases.  This  is  becuase  the  target  signal  dominates  in  these 
cases.  On  the  other  hand  for  the  2000m  and  5000m  cases  the  optimal 
filter  configuration  shows  a  slight  increase  in  tracking  performance, 
than  the  stab  error  not  fed  forward  configuration.  Such  a  conclu¬ 
sion  is  reached  by  noting  that  the  random  "spikes"  of  the  optimal 
filter  configuration  are  less  in  magnitude  than  for  the  stab  error 
not  feed  forward  configuration.  The  reason  for  this  increase  in 
performance  is  that  BMT  plays  a  larger  role  in  the  beam  error  for 

these  cases  and  the  stab  error  is  fed  forward  as  part  of  the  scheme 
to  minimize  the  effects  of  BMD  on  the  beam  error. 


i  P  r~i  MntiV  r'~Yi 'mTiii" 
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In  order  to  determine  how  the  power  of  the  error  signal  is 
distributed  with  respect  to  frequency  the  power  spectrum  density 
of  each  of  the  signals  shown  in  Figures  9-13  were  generated 
using  an  HP-85  desktop  computer  and  the  standard  Waveform  Analysis 
Pac  purchased  from  HP  [Ref,  5].  The  PSD  of  each  error  signal  is 
shown  in  Figures  14-18.  The  frequency  range  in  each  case  is 
0-50  Hs.  It  should  be  noted  that  the  ordinate  axis  is  in 
declbels/Hs  but  has  been  normalised  as  indicated. 

For  the  Figures  14-15,  it  is  seen  that  most  of  the  power 
occurs  in  the  low  frequency  range.  This  is  caused  by  the  fact 
that  the  tracker  error  is  not  fed  to  the  AAS  in  any  of  these 
cases,  and,  hence,  the  capability  of  the  system  to  track  a 
moving  target  is  limited. 

On  the  other  hand  in  Figures  16-18,  it  is  seen  that  the  error 
power  is  more  evenly  spread  across  the  spectrum  and  tends  to  roll¬ 
off  on  the  high  frequency  end.  This  was  expected  since  it  is 
known  that  the  PSD’s  of  the  BMD  tracker  error  and  target  signal 
roll-off  at  the  higher  frequencies. 

In  closely  scrutinizing  Figures  16-18,  only  in  the  500m  cases 
is  any  resonance  effects  detected.  A  clear  resonance  is  seen  to 
occur  at  about  40  Hz  (about  204.8  on  the  normalized  frequency 
scale  of  Figures  16-18) .  This  is  attributed  to  the  target  having 
enough  motion  that  natural  resonances  of  the  stabilized  platform 
on  which  the  telescope  is  mounted  are  excited  and,  hence  contributed 
significantly  to  the  beam  error. 
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FREO.  INTERVAL*  .1953125  ©  FREQ.  INTERVAL®  1953125 

MAX  VALUE  =  8055.68817026  -  MAX  VALUE  *  410.282529708 

(c)  2000m  (d)  joooa 

Figure  14.  PSD  of  Beam  Error  for  No  Feed  Forward  Configuration 
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Figure  18.  PSD  of  Beam  Error  for  Optimal  Filter  Configuration 


The  Figures  also  verify  that  the  configurations  where  the 
optimal  filters  (Figure  16  and  18)  were  used,  produced  less  beam 
error  than  the  unity  filter  configuration  (Figure  17),  One 
measure  of  this  is  gained  by  comparing  the  maximum  values  of  the 
corresponding  cases.  By  similar  comparisons  it  is  seen  that  the 
optimal  filter  configuration  with  the  stab  error  fed  to  the  IMC 
and  to  the  AAS  (Figure  18)  is  slightly  superior  to  the  stab  error 
not  feed  forward  case  (Figure  16). 

Figures  9-18  show  that  much  can  be  gained  from  the  feed 
forward  scheme.  A  main  question  is  whether  the  optimal  filter 
schemes  produce  enough  improvement  in  order  to  justify  the 
additional  electronics.  As  was  pointed  out  in  [3]  and  earlidr 
in  this  report,  the  logical  choice  for  P(s)  is  unity  if  an 
optimal  P(s)  is  not  used.  In  order  to  provide  a  clear  comparison 
of  the  unity  filter  and  optimal  filter  configurations  for  each 

case  of  the  flyby  type  scenario  the  PSD's  of  each  were  numerically 

* 

integrated  backward  starting  at  50  Ha.  The  corresponding  cases 

are  shown  in  Figures  19-20,  where  the  ordinate  values  have  been 
-12 

scaled  by  10  .  In  interpretating  each  plot  an  ordinate  value 

gives  the  total  amount  of  power  in  all  frequencies  greater  than 
to  the  corresponding  frequency  (up  to  50  Hz.),  and  the  difference  in 
any  two  ordinate  values  gives  the  power  for  the  corresponding 
frequency  range. 


*Very  little  power  is  contained  in  the  spectrum  for  frequencies 
greater  than  50  Hz. 
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Figure  19,  Backward  Sum  of  Unity  Feed  Forward/Optimal 
Configurations 


Figure  19.  (Continued)  Backward  Sum  of  Unity  Feed  Forward/ 

Optimal  Configurations 


If  a  point  by  point  comparison  of  corresponding  cases  is 
done,  it  is  easily  seen  that  the  optimal  filter  configuration 
cases  are  in  every  instance  below  the  unity  filter  configuration 
cases.  This  verifies  that  the  optimal  filters  do  reduce  the 
error  signals. 

It  should  also  be  noted  that  the  backward  sum  plots 
definitely  verify  that  a  resonance  mode  of  the  system  was  excited 
in  the  500m  case.  This  is  seen  by  a  sharp  rise  in  the  backward 
sum  curves  at  about  40  Hz  for  the  500m  cases.  In  addition  the 
backward  sum  curves  also  show  the  possibility  of  another  resonance 
being  excited  at  about  19  Hz  for  the  500m  case.  Apparently  in 
none  of  the  other  cases  were  these  resonances  excited  enough  to 
contribute  significantly  to  the  beam  error. 

The  backward  sum  plots  can  also  be  used  to  obtain  RMS  errors 
for  the  beam  angle.  Since  these  are  a  backward  integration  of 
two-sided  PSD's  then 

RMS  ERROR  *  ^  2  x  ORDINATE  AXIS  INTERCEPT 

For  the  cases  at  hand  the  RMS  errors  are  shown  in  Table  II.  Again 
this  shows  that  the  optimal  filter  can  definitely  improve  tracking 
performance.  Furthermore,  it  should  be  noted  that  the  RMS  values 
in  Table  II  are  comparable  to  those  obtained  in  [3].  Actually  the 
RMS  values  are  slightly  higher  than  those  obtained  in  [3].  This  is 
expected  since  the  sampling  effects  of  the  tracker  are  included  in 
this  study  and  models  of  tracker  noise  and  BMD  that  covered  a  wider 
spectrum  are  used. 
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Yo 

UNITY  FILTER 
RMS  ERROR 
prad 

OPTIMAL  FILTER 

RMS  ERROR 
prad 

500m 

8.94 

7.62 

1000m 

5.10 

4.69 

2000m 

4.96 

4.00 

5000m 

4.96 

3.87 

TABLE  II:  COMPARISON  OF  RMS  ERRORS  FOR 
UNITY  FILTER  vs  OPTIMAL 
FILTER  CONFIGURATIONS 
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V.  CONCLUSION  AND  RECOMMENDATIONS 


A  feed  forward  scheme  for  controlling  the  LOS  of  a  high  energy 
laser  beam  was  devised  for  reducing  errors  of  the  beam  angle  due  to 
BMD  and  for  tracking  targets  with  significant  motion.  As  was 
shown  in  [1]  and  [2]  the  capability  of  the  system  to  reject  BMD 
was  greatly  enhanced  by  this  scheme,  but  the  effects  on  beam  error 
due  to  tracker  noise  was  amplified.  A  third  study  [3]  was  initiated 
to  design  a  compromising  filter  that  would  provide  the  capabilities 
to  satisfactorily  reject  BMD  and  tracker  noise  signals  and  to 
track  targets  with  significant  motion.  This  was  accomplished  by 
designing  a  filter  to  minimize  the  RMS  error  of  the  beam  angle. 

However,  the  study  in  (3]  used  RMS  values  and  simplifying 
assumptions  in  regard  to  the  system,  e.g.,  the  sampling  effects 
of  the  tracker  were  neglected.  Since  the  results  produced  by  [3] 
appeared  very  promising  it  was  decided  that  follow-on  studies 
should  be  initiated.  In  particular,  it  was  decided  that  time 
domain  studies  should  be  conducted  to  verify  the  previous 
frequency  domain  studies. 

These  time  studies  are  summarized  in  this  report.  In 
accomplishing  them  some  of  the  simplifying  assumptions  made  in 
[3]  were  removed,  e.g.,  the  sampling  effects  of  the  tracker  were 
included.  The  results  contained  in  this  report  closely  agree  with 
those  contained  in  [3].  As  a  general  rule  the  RMS  beam  angle 
errors  are  slightly  greater,  but  this  was  expected  since  more 


exact  models  and  signals  were  used.  Ad  a  result  of  the  studies  In 

[1],  [2],  [3],  and  this  repo  t  the  feed  forward  scheme  appears  to 
be  very  promising  for  beam  control. 

There  are  several  follow-on  studies  that  should  be  conducted. 
Among  these  are 

(1)  The  scheme  has  been  tested  with  an  ideal  flyby  type 
scenario.  It  should  be  tested  with  typical  scenarios. 

(2)  The  study  in  [3]  produces  different  filters  for  e.ich  case 
of  the  flyby  type  scenario.  These  designs  should  be 
repeated  for  other  scenarios,  and  it  should  be  determined 
if  a  pattern  exists  for  determining  the.  filter  coefficients 
for  a  general  type  scenario.  It  might  very  well  turn  out 
that  an  adaptive  filter  could  be  used,  i.e.  based  upon 
target  motion  the  system  should  automatically  update 

the  coefficients  of  the  filter. 

(3)  More  accurate  models  of  the  system  signals  and  components 
should  be  used  in  the  design  of  the  optimal  filters. 

In  accomplishing  any  or  all  of  the  above  suggested  studies 
the  code  developed  in  [3]  (or  a  modified  version)  could  be  used. 

Of  course  each  of  the  suggested  frequency  domain  studies  should  be 
paralleled  with  time  domain  studies  to  verify  that  the  system  is 
operating  within  its  designed  dynamical  range. 
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APPENDIX 


The  following  is  a  listing  of  the  APT/AAS  MIMIC  Simulation 
with  corresponding  data  for  the  500m  case.  Definitions  of  several 
of  the  pertinent  constants  and  parameters  used  in  the  code  are  aa 
follows: 

SPD  -  Sampling  period  of  tracker 
SPD2  -  Sampling  period  used  for  generating  BMD 
V  -  Target  Velocity  (m/s) 

YO  -  Minimum  Approach  distance  s(m) 

TDL  -  Tracker  Delay  (s) 

LCV12,  LCV13,  &  LCV14  -  Logical  Control  Variables  for  Determining 
the  Control  Configuration  (1.  -  True,  0.  -  False) 

All,  A12,  Bll,  and  B12  -  Numerator  Coefficients  of  Optimal  Filter 
A21,  a22,  B21,  and  B22  -  Denominator  Coefficients  of  Optimal  Filter 
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onnnornnonnonnnnonoo  oo  ooo 


SIMULATION  OF  THE  BEAM  CONTROL  SYSTEM  USING 
BOTH  THE  APT  AND  AAS  TO  FORM  VARIOUS  FEEDFORWARD 

SCHEMES 


CON<TZ,SPD,TBMD,SPD2> 

CON(Ct ,Tt ,C2,T2,C3,C4> 

CON < VAR 1 , HUH t , V AR 2 , NUM2 > 

CON ( SIGNAL , TR  K  ER R , BMD  > 

CON ( DT , DTMAX  ,DTM IN, PERIOD, TMAX ) 
CON<V,YQ  > 

CON<A3,A2,At,AO> 

C0N<P1,P2,P3,P4,P5,PA> 

CON<P7,Pe,P*,PiO,Ptl,Pi2> 

CON<TDL,TDX,Xl> 

C0N<X4,X66,X6> 

CON ( LCV1 2 , LCVt  3 , LCV 14) 
PAR<Att,A12,Blt ,B12> 
PAR<A2t,A22,B21,B22> 


IF  LCVt 2=*TRUE , THEN  THE  TRACK  ERROR  IS  FED  FORWARD. 

IF  LCVt 2-FALSE, THEN  THE  TRACK  ERROR  IS  NOT  FED  FORWARD. 

IF  LCVt 3-TRUE, THEN  THE  STAB  ERROR  IS  FED  TO  THE  AAS. 

IF  LCVt 3-FALSE, THEN  THE  STAB  ERROR  IS  NOT  FED  TO  THE  AAS. 

IF  LCVt 2  AND  LCVt 3  ARE  TRUE, THE  COMPLETE  FEEDFORWARD 
SCHEME  IS  IMPLEMENTED. 

IF  LCV12  AND  LCVt 3  ARE  FALSE, THE  BEAM  IS  CQNTROLED  BY  THE  APT  0 


IF  LCVt 4-TRUE, THEN  THE  UNITY  TRACK  ERROR  FEEDFORWARD 
FILTER  IS  IMPLEMENTED. 


GENERATION  OF  TARGET  SIGNAL 


TX 

T-6 . 

LCVt 

FSWCT, TRUE, TRUE, FALSE) 

LCV2 

FSW(T-4. , TRUE , TRUE jFALSE > 

LCV3 

NOT1LCV2) 

LCV4 

FSW<T-8. , TRUE, TRUE, FALSE) 

LCVS 

AND  <  LCV3 , LCV4  > 

LCV6 

NOT<LCV4> 

LCV7 

FSW ( T-TMAX , FALSE , FALSE , TRUE ) 

LCV13 

Ptt 

t .  0 

LCVt 

com 

324 . #A3-72 . «A2+t  8 . *At -6 . »A0 

LCVt 

CONS 

-4 . *A3+8 . *A2/3 . -2 . «At  +2 . *A0 

LCVt 

COM3 

2, *V/Y0 

LCV1  XP2  v#a. 

LCV1  MIT  -2. 

LCVt  XH2  V#MIT 

LCVt  CONS  (((♦A3/4,*MXT*-A2/3.  )#MIT+At/2.  )*MITt>AQ  )#MXT-CONl 
LCVt  L0N6  ATN(XM2/Y0) 

LCV3  X  V»TX 

LCV5  THETA  ATN(X/YO)-CON6+CON5 

LCV2  PARTI  (((♦A3/4.Q*TX+A2/3.Q)*TX+At/2,Q)»TX*AO)#TX-CON1 

LCVt  C0N4  ATN( XP2/Y0 ) -CQN&+C0N5-C0N2 

LCV6  PART2  ( <  (-A3/4 . 0*TX+A2/3, 0)*TX-Al/2. 0)»TX+AQ)»TX+CON4 

SIG1  LSW(LCV2,PARTt ,0. ) 

SIC2  LSW(LCV5, THETA, 0. > 

SIG3  LSW(LCV6,PART2,0. ) 

SIGNAL  SIGH-SIG2+SIG3 

GENERATION  OP  BHD 

LCV1 0  FSW  <  T-TBMD , FALSE , TR  UE , TR  UE ) 

LCVtO  BHD I  C2«RNG(0. ,VAR2,NUM2) 

LCV10  TBHD  TBHD+SPD2 

BHD2  I NT (BHDI >0.0) 

BHD3  C4*FTR (BMDI , T2) 

BHD  C3* ( BMD2+BMD3 > 

IBHD  INT (BHD >  0 , > 

ISBHD  INT  <  BMD*BMD ,  0  .  ) 

TRACKER  ERROR  SIGNAL 

ICV9  FSW ( T-TZ , FALSE , TRUE >  TRUE ) 

LCV9  TRKERR  C1*RNG( 0 . , VARt ,NUMt > 

LCV9  TZ  TZ+SPD 

INTRER  INT (TRKERR ,0 . ) 

INSQTE  INT (TRKERR*TRKERR , 0 . > 

SIG  ZQH( SIGNAL, SPD) 

LCVt  2  X66  Z0H(X6,SPD) 

GANG1  ZOH(UANG,SPD> 

OUTPUT  OF  HOLD  DEVICE 

XI t  SIG+TRKERR-Pt 1*X66-GANG1 

OUTPUT  OF  DELAY 

LCV1 1  FSW(T-TDX, FALSE, TRUE, TRUE) 

LCVt  l  TDX  TDX+SPD 

LCVtt  XI  XI t 

OUTPUT  OF  TRACK  LOOP  COHPENSATOR 

X2  33 . 02#X 1+350 .0*1 NT (XI, 0.0) 

STABILIZE  PLATFORM  ERROR  SIGNAL 
ESTAB  INT ( X2-X3 ,0.0) 

PLATFORM  SIGNAL 

X3  BHD-tWl 

PLATFORM  DYNAMICS 

W1  INT ( U2 , 0 . 0  ) 

W2  INT <U3, 0 . 0  ) 

W3  INT(Pl#W3+P2*W2+P3«Wt+P10#Z1  ,0.0) 

Z1  P4*ESTAB+P5*INT (ESTAB, 0.0) 

GIHBAL  ANGLE 

GANG  INT <X3, 0.0) 
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OPTIMAL  FILTER  TRANSFER  FUNCTION 
LCV16  N0T<LCV14> 

LCV15  AND<L.CV12,LCV16> 

LCV1 4  X4  XI 

LCV15  X41  X1*A11/B1H-FTR<<A12/B12~A11/B11 >#Xt ,B11/B12> 

LCV15  X4  X4t  #A21/B21+FTR ( <A22/B2°-A21/B21 >#X41 , B21/B22) 

100  HZ  STAB  FILTER 
LCV13  X6  FTR  <ESTAB ,F 1 2 ) 

AUTO  ALIGNMENT  SYSTEM  ERROR 

EAAS  X4+P  1 1 #X6-X5 

AUTO-ALIGNMENT  SYSTEM  DYNAMICS 

X51  P6#£AAS+P7#INT (EAAS , 0 . 0 ) 

X5  INT  <  X51 , 0 . 0 ) 

BEAM  ANGLE 

BANG  X5+GANG 

ERROR  BETWEEN  BEAM  ANGLE  AND  TARGET  SIGNAL 
ERROR  '  SIGNAL-BANG 
ROOT  MEAN  SQUARE  VALUE  OF  ERROR 

MSE  I NT<ERROR TERROR /PERIOD, 0 . 0) 

HDR <  TIME , SIGNAL ,GIMBAL , BEAM , ERROR ,MSE> 

HDR (  ,  , ANGLE, ANGLE) 

HDR 

HDR  <  ,I.TRK,IS.TRK,I. BMD , IS . BMD) 

HDR 

OUT ( T , SIGNAL , GANG , BANG , ERROR , MSE ) 

OUT  <  , INTRER , INSQTE , IBMD, ISBMD) 

OUT 

FIN(T,TMAX> 

END 

0.  ,01666663666  0.  6.28319  E-03 

2.192  E-06  2.652  E-03  -1.0  .005  9.403  E-04  .1061 

2.0  3975319.  2.0  1234567. 

1.0  0.0  0.0 

.01  3.0  E-03  1.0E  E-05  8.0  8.0 

457.2  5000.0 

-2 . 586056E-3-3 . 067474E-2-8 . 880281E-21 . 288565E-02 


-279.624 

-369514.11 

-78369905.0  265.3 

10000.0 

1987. 

1.248  E 

06.0045475 

-352692.0  78369905. 

0,0 

1 . 592E 

.003 

.0  03 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 .3745 

7 . 9788 

1,6585  10.0 

1.1542 

100.09 

2.6441  80.0 

ENDATA 


STOP 
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OPTIMAL  FILTER  TRANSFER  FUNCTION 
LCV16  NOT (LCV14> 

LCV15  AND(L.CV12,LCV16> 

LCM14  X4  XI 

LCV15  X41  XI *A1 1 /B1 1+FTR  < (A12/B12-A11/B11 )#Xt , B1 1/B12) 

LCVtS  X4  X41*A21/B2H-FTR((A22/B22-A21/B21)*X41>B21/B22> 

100  HZ  STAB  FILTER 

LCV13  X6  FTR<ESTAB,P12> 

AUTO  ALIGNMENT  SYSTEM  ERROR 

EAAS  X4«-P11*X6-X5 

AUTO-ALIGNMENT  SYSTEM  DYNAMICS 

X51  P6#EAAS+P7#INT  <EAAS  >0.0) 

X5  INT (X51 , 0 . 0 ) 

BEAM  ANGLE 

BANG  X5+GANG 

ERROR  BETWEEN  BEAM  ANGLE  AND  TARGET  SIGNAL 
ERROR  SIGNAL-BANG 

ROOT  MEAN  SQUARE  VALUE  OF  ERROR 

MSE  INT (ERROR TERROR /PER IOD,  0 . 0 ) 

HDR (TIME, SIGNAL , GIMBAL , BEAM , ERROR ,MSE  > 

HDR (  ,  , ANGLE, ANGLE) 

HDR 

HDR(  ,I.TRK,IS.TRK,I. BMD , IS . BMD) 

HDR 

OUT ( T , SI GNAL . GANG , BANG , ER  ROR „ MSE ) 

OUT (  , INTRER  INSQTE , IBMD , ISBMD) 

OUT 

FIN<  T ,TMAX) 

END 

0.  . 01666666666  0,  6.28319  E-Q3 

2.192  E-06  2.652  E-03  -1.0  .005  9.403  E-04  .1061 

2.0  3975319.  2.0  1234567. 

1.0  0.0  GO 

.01  3.0  E-03  1.0E  E-05  8.0  8.0 

457.2  5000.0 

-2 . 586056E-3-3 . 067474E-2-8 . 880281 E-21 . 288565E-02 

-279.624  -369514.11  -78369905.0  265.3  10000.0  1987. 

1.248  E  06.0045475  -352692.0  78369905.  0.0  1.592E-03 

.  003  .003  0.0 

0.0  0.0  0.0 

0.0  C.O  0,0 

1.3745  7. 9788  1.6585  10.0 

1.1542  100.09  2.6441  80.0 

LNDATA 

STOP 

EOF i 149  SCAN; 44 
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